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Abstract 

The results of background measurements with the second version of the BOREX- 
INO Counting Test Facihty (CTF-II), installed in the Gran Sasso Underground 
Laboratory, were used to obtain limits on the instability of nucleons, bounded in 
nuclei, for decays into invisible channels {inv): disappearance, decays to neutrinos, 
etc. The approach consisted of a search for decays of unstable nuclides resulting 
from A'" and NN decays of parent ^^C, ^^C and ^^O nuclei in the liquid scintilla- 
tor and the water shield of the CTF. Due to the extremely low background and 
the large mass (4.2 ton) of the CTF detector, the most stringent (or competitive) 
up-to-date experimental bounds have been established: T(n inv) > 1.8 • 10^^ y, 
t{p inv) > 1.1 • 10^^ y, T{nn inv) > 4.9 • 10^^ y and T{pp — > inv) > 5.0 • 10^^ 
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y, all at 90% C.L. 



Key words: proton decay, baryon number conservation, scintillation detector 
PACS: 11.30, 11.30.F, 12.60, 29.40.M 



1 Introduction 

The baryon (B) and lepton (L) numbers are considered to be conserved in the 
Standard Model (SM) ^ . However, no symmetry principle underlies these laws, 
such as, e.g. gauge invariance, which guarantees conservation of the electric 
charge. Many extensions of the SM include B and L violating interactions, 
predicting the decay of protons and neutrons bounded in nuclei. Various decay 
mechanisms with AB—1, 2 and A{B-L)—0, 2 have been discussed in the 
literature intensively [2,3]. A novel baryon number violating process, in which 
two neutrons in a nucleus disappear, emitting a bulk majoron nn x, was 
proposed recently [4]; the expected mean lifetime was estimated to be ~10"^^^^^ 
y. Additional possibilities for the nucleon (iV) decays arc related to theories 
which describe our world as a brane inside higher-dimensional space [5,6]. 
Particles, initially confined to the brane, may escape to extra dimensions, thus 
disappearing for the normal observer; the characteristic proton mean lifetime 
was calculated to be r(p)=9.2-10^'' y [7]. Observation of the disappearance of 
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e , N, NN would be a manifestation of the existence of such extra dimensions 
[6]. 

No evidence for nucleon instabihty has been found to date. Experimental 
searches [8] with the 1MB, Prejus, (Super) Kamiokande and other detectors 
have been devoted mainly to nucleon decays into strongly or electromagnet- 
ically interacting particles, where lower limits on the nucleon mean lifetime 
of 10'^°^^^ y were obtained [9]. At the same time, for modes where N or A^A^ 
pairs disappear or they decay to some weakly interacting particles (neutrinos, 
majorons, etc.), the experimental bounds are a few orders of magnitude lower. 
Different methods were applied to set limits for such decays ^ (see table 1 for 
summary) : 

(1) Using the hmit on the branching ratio of spontaneous fission of '^^'^Th 
under the assumption that p or n decay in ^^^Th will destroy the nucleus 
[10]. The bound on the mean lifetime obtained in this way can be considered 
independent of the a p or n decay mode, since the ^^^Th nucleus can be 
destroyed either by strong or electromagnetic interactions of daughter particles 
with the nucleus or, in the case of N disappearance, by subsequent nuclear 
deexcitation process; 

(2) Search for a free n created after p decay or disappearance in the deuterium 
nucleus {d=pn) in a liquid scintillator enriched in deuterium [11] or in a volume 
of D2O [12,16,17]; 

(3) Geochemical [13] or radiochemical [14] search for daughter nuclides which 
have appeared after A^ decays in the mother nuclei (valid for decays into 
invisible channels); 

(4) Search for prompt 7 quanta emitted by a nucleus in a de-excitation process 
after A^ decays within the inner nuclear shell [18] (valid for invisible channels); 

(5) Considering the Earth as a target with nucleons which decay by emitting 
electron or muon neutrinos; the z/g, t'^ can be detected by a large underground 
detector [19,20] (valid for decay into neutrinos with specific fiavors); 

(6) Search for bremsstrahlung 7 quanta emitted due to a sudden disappearance 
of the neutron magnetic moment [21] (limits depend on the number of emitted 
neutrinos) ; 

(7) Study of radioactive decay of daughters (time- resolved from prompt prod- 



We are using the following classification of decay channels: decay to invisible 
channel (inv) means disappearance or decay to weakly interacting particles (one 
or few neutrinos of any flavors, majorons, etc.). Channel to anything, mentioned 
below, includes also decays to inv. 
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Table 1 

Lower limits on the mean lifetime for decay of nucleons, bounded in nuclei, into 
invisible channels established in various approaches (see footnote 7 for channels 
classification). 



N or NN decay 



r limit, y 
and C.L. 



Year, reference and short explanation 



P 



n 



nn 



pp 
pn 



anything 


1.2-1023 




1958 




3.0-1023 




1970 




4.0-1023 


95% 


2001 


inv 


7.4-1024 




1977 




1.1-1026 




1978 




1.9-1024 


90% 


2000 




~1028 




2002 




3.5-1028 


90% 


2003 




1.1-1026 


90% 


2003 


anything 


1.8-1023 




1958 


inv 


8.6-1024 




1977 




1.1-1026 




1978 




4.9-1026 


90% 


1993 




1.8-1025 


90% 


2003 




5.0-1026 


90% 


1979 




1.2-1026 


90% 


1991 




3.0-1025 


90% 


1991 


ViViVi 


2.3-102^ 


90% 


1997 


ViViViViVi 


1.7-1027 


90% 


1997 




6.0-1024 


90% 


1991 




1.2-1025 


90% 


1991 


inv 


1.2-1025 


90% 


2000 




4.2-1025 


90% 


2002 




4.9-1025 


90% 


2003 


inv 


5.5-1023 


90% 


2000 




5.0-1025 


90% 


2003 


inv 


2.1-1025 


90% 


2002 



10] limit on spontaneous fission 

11] search for free n in liquid scintillator 

enriched in deuterium [d n+?) 
12] free n in D2O volume 
13] geochem. search for ^30Te^...^^29xe 
14] radiochem. search for ^'^K^...^'^^ Av 
15] search for ^2^1 decay in ■'^29xe detector 
16] free n in D2O volume 
17] free n in D2O volume 
" ] search for ^'^B decay in the CTF detector 
10] limit on ^'^^Tli spontaneous fission 
13] geochem. search for ^30^6^...— ^^^^Xe 
14] radiochem. search for ^^K^...—^"^"^ At 
18] search for 7 with £"^,=19-50 McV emitted in 

^50 deexcitation in Kamiokande detector 
] search for ^^C decay in the CTF 
19] massive liquid scint. detector fired by U/j, 

in result of n decays in the whole Earth '''^ 
20] Frcjus iron detector fired by i^^ 
20] Frejus iron detector fired by Vg 
21] search for bremsstr. 7 with £'-y>100 MeV 

emitted due to sudden disapp. of n magn. 

moment (from Kamiokande data) ^ 
21] the same approach 
20] Frejus iron detector fired by ^ 
20] Frejus iron detector fired by ^ 
15] search for ■'^2'''xe decay in ^29xe detector 
22] radiochem. search for 39K— >...— >37^j. h 
" ] search for -'^''C and "'^40 decay in the CTF 
15] search for ^'^'^Te decay in ^29Xe detector 
" ] search for ^^Be decay in the CTF detector 



22] radiochem. search for 39K- 



.37 Ar 



" This work 

^ The result of [19] was reestimated in [20] to be more than one order of magnitude lower 
^ The limit is also valid for p — > fpyi^u^ decay 
The limit is also valid for p Ue^e^e decay 

^ i = e, /i, r 

^ The limit is also valid for pn and pp decays into VfjV^ 
^ The limit is also valid for pn and pp decays into u^pe 
^ On the base of the data of ref. [14] 
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ucts), created as a result of N or NN decays of the mother nuclei, incorporated 
into a low-background detector (valid for decay into invisible channels). This 
method was first exploited by the DAMA group with a liquid Xe detector [15]. 

In the present paper we use the same approach to search for N and NN 
instability with the Counting Test Facility, a 4.2 ton prototype of the multiton 
BOREXINO detector for low energy solar neutrino spectroscopy [23]. The 
preliminary results were presented in [24]. 



2 Experimental set-up and measurements 

2.1 Technical information about CTF and BOREXINO 

BOREXINO, a real-time 300 ton detector for low-energy neutrino spectroscopy, 
is nearing completion in the Gran Sasso Underground Laboratory (see [23] and 
refs. therein). The main goal of the detector is the measurement of the ^Be 
solar neutrino flux via v — e scattering in an ultra-pure liquid scintillator, 
while several other basic questions in astro- and particle physics will also be 
addressed. 

The Counting Test Facihty (CTF), installed in the Gran Sasso Underground 
Laboratory, is a prototype of the BOREXINO detector. Detailed reports on 
the CTF results have been published [23,25,26], and only the main character- 
istics of the set-up are outlined here. 

The CTF consists of an external cyhndrical water tank (011x10 m; RilOOO 
t of water) serving as passive shielding for 4.2 m^ of hquid scintillator (LS) 

contained in an inner spherical vessel of 02.0 m. High purity water with a 
radio-purity of ^lO^^^ g/g (U, Th), ^IQ-^"^ g/g (K) and <2 /xBq// for 222Rn 
is used for the shielding. The LS was purified to the level of ~10~^^ g/g in 
U/Th contamination. 

We analyze here the data following the upgrade of the CTF (CTF-II). The 
liquid scintillator used at this stage was a phenylxylylethane (PXE, CieHig) 
with p-diphenylbenzene (para-terphenyl) as a primary wavelength shifter at 
a concentration ol2 g/l along with a secondary wavelength shifter l,4-bis-(2- 
methylstyrol)-benzene (bis-MSB) at 50 m.g/1. The density of the scintillator is 
0.996 kg//. The scintillator principal de-excitation time is less than 5 ns which 
provides a good position reconstruction. In the CTF-II an additional nylon 
screen between the scintillator vessel and PMTs (against radon penetration) 
and a muon veto system were installed. 
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The scintillation light is collected with 100 phototubes (PMT) fixed to a 7 
m diameter support structure inside the water tank. The PMTs are fitted 
with light concentrators which provide a total of 21% optical coverage. The 
number of photoelectrons measured experimentally is 3.54 per PMT for 1 MeV 
electrons at the detector's center. 

For each event the charge and timing of hit PMTs are recorded. Each channel is 
supported by an auxiliary channel, the so-called second group of electronics, 
used to record all events coming within a time window of 8.2 ms after the 
trigger, which allows tagging of fast time-correlated events with a decrease of 
the overall dead time of the detector. For longer delays, the computer clock 
is used providing the accuracy of f^iO.l s. Event parameters measured in the 
CTF include the total charge collected by the PMTs during 0-500 ns, used to 
determine an event's energy; the charge in the "tail" of the pulse (48-548 ns) 
which is used for pulse shape discrimination; PMT timing, used to reconstruct 
the event's position; and the time elapsed between sequential events, used to 
tag time-correlated events. 



2.2 Detector calibration 

The energy of an event in the CTF detector is defined using the total collected 
charge from all PMT's. In a simple approach the energy is supposed to be lin- 
ear with respect to the total collected charge. The coefficient finking the event 
energy and the total coHected charge is called light yield (or photoelectron 
yield). The light yield for electrons can be considered linear with respect to 
its energy only for energies above 1 MeV. At low energies the phenomenon of 
"ionization quenching" violates the linear dependence of the light yield ver- 
sus energy [27]. The deviations from the linear law can be taken into account 
by the ionization deficit function f{kB,E), where ks is an empirical Birks' 
constant. For the calculations of the ionitazion quenching effect for PXE scin- 
tillator we used the KB program from the CPC library [28]. The "ionization 
quenching" effect leads to a shift in the position of the full energy peak for 
gammas on the energy scale calibrated using electrons. In fact, the position 
of the 1461 keV ''^K gamma in the CTF-II data corresponds to 1360 keV of 
energy deposited for an electron. 

The detector energy and spatial resolution were studied with radioactive sources 
placed at different positions inside the active volume of the CTF. A typical 
spatial la resolution is 10 cm at 1 MeV. The studies showed also that the total 
charge response of the CTF detector can be approximated by a Gaussian. For 
energies E>1 MeV (which are of interest here), the relative resolution can be 

expressed as as/E = ^3.8 • lO-^/E + 2.3 • 10-3 {E is in MeV) [29] for events 
uniformly distributed over the detector's volume. 
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The energy dependence on the coUected charge becomes non-hnear for the 
energies ^ 4.5 MeV in the first group of electronics because of the saturation 
of the ADCs used. In this region we are using only the fact of observing or not- 
observing candidate events, hence the mentioned nonlinearity doesn't influence 
the result of the analysis. 

Further details on the energy and spatial resolutions of the detector and ion- 
ization quenching for electrons, 7 quanta and a particles can be found in 
[29,30]. 

2.3 Muon veto 

The upgrade of the CTF was equipped with a carefully designed muon veto 
system. It consists of 2 rings of 8 PMTs each, installed at the bottom of the 
tank. The radii of the rings are 2.4 and 4.8 m. Muon veto PMTs are looking 
upward and have no light concentrators. The muon veto system was optimized 
in order to have a negligible probability of registering the scintillation events 
in the so-called "neutrino energy window" (250-800 kcV). The behaviour of 
the muon veto at the higher energies has been specially studied for the present 
work. Experimental measurements with a radioactive source (chain of ^^^Ra) 
[31] gave, for the probability ri{E) of identification of an event with energy E 
in the LS by the muon veto, the value of (1±0.2)% in the 1.8-2.0 MeV region 
(see fig. 3 later). The energy dependence of 1]{E) was also calculated by a 
ray-tracing Monte Carlo method accounting for specific features of the light 
propagation in the CTF which are detailed in [26]. The calculated function 
was adjusted to reproduce correctly the experimental measurements with the 
^^^Ra source. 

2.4 Data selection 

As it will be shown below, the candidate events, relevant for our studies, have 
to satisfy the following criteria: (1) the event should occur in the active volume 
of the detector and must not be accompanied by the muon veto tag; (2) it 
should be single (not followed by a time-correlated event); (3) its pulse shape 
must correspond to that of events caused by 7 or /3 particles. 

The selection and treatment of data (spatial cuts, analysis of an event's pulse 
shape to distinguish between electrons and a particles, suppression of external 
background by the muon veto system, etc.) is similar to that in ref. [30]. 

The experimental energy spectra in CTF-II, accumulated during 29.1 days of 
measurements, are shown in fig. 1. The spectrum without any cuts (spectrum 
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1) is presented on the top. The second spectrum is obtained by applying 
the muon cut. which suppressed the background rate by up to two orders of 
magnitude, depending on the energy region. 




Fig. 1. Background energy spectra of the 4.2 ton BOREXINO CTF-II detector 
measured during 29.1 days. From top to bottom: (1) spectrum without any cuts; 
(2) with muon veto; (3) only events inside the radius i?<100 cm with additional 
a/ (3 discrimination applied to eliminate any contribution from a particles; (4) pairs 
of correlated events (with time interval At<8.2 ms between signals) are removed. 
In the inset, the simulated response function for external ^'^K gammas is shown 
together with the experimental data. 

On the next stage of the data selection we applied a cut on the reconstructed 
radius. In the energy region 1-2 MeV we used R <100 cm cut aiming to 
remove the surface background events (mainly due to the decays outside 
the inner vessel) and leave the events uniformly distributed over the detector 
volume. The efficiency of the cut has been studied with MC simulation and 
lays in the range of tR = 0.76 — 0.80 in the energy region 1-2 MeV. Additional 
a/f3 discrimination [25] was applied to eliminate contribution from a particles 
(spectrum 3 in fig. 1). 
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The time-correlated events (that occurred in the time window At<8.2 ms) 
were also removed (spectrum 4). The peak at 1.46 MeV, present in all spec- 
tra, is due to decays outside the scintillator, mainly in the ropes supporting 
the nylon sphere. The peak-hke structure at ~6.2 MeV is caused by satura- 
tion of the electronics by high-energy events. The lower spectrum of fig. 1 
presents all the candidate scintillation events in the search for the decay of 
radioactive nuclides created in the active volume of the CTF after the nucleon 
disappearance in parent nuclei. 



3 Data analysis and results 



3.1 Theoretical considerations 



The decay characteristics of the daughter nuclides, resulting from and NN 
decays in parent nuclei - ^^C, ^^C and ^^O - contained in the sensitive volume 
of the CTF liquid scintillator or in the water shield, are listed in table 2. 

Table 2 

Initial nuclei in the Counting Test Facility, their abundance 8 [32], processes with 
A5 = 1,2 and characteristics of daughter nuclides [33]. 



Initial De- Daughter nucleus, half-life, modes of decay and energy release 

nucleus cay 



12p 

5=98.93% 



5=1.07% 



16q 



n 
P 

nn 

pn 
pp 

n 
P 

nn 
pn 
pp 

n 



'IB 

10 p 

10b 



lip 

'IB 



15 



o 



5=99.757% p 



nn 
pn 
pp 



14 Q 

14 p 



ri/2=20.38 m 
stable 
Ti/2=19.2 s 

stable 

ri/2=1.6-106 y 



stable 

Ti/2=20.4 ms 
ri/2=20.38 m 
stable 
Ti/2=13.8 s 



ri/2=122 s 
stable 

ri/2=70.60 s 
stable 

Ti/2=5730 y 



P+ (99.76%), EC(0.24%); Q=1.982 MeV 
/?+; Q=3.651 MeV 
p-; Q=0.556 MeV 

Q=13.370 MeV 
P+ (99.76%), EC (0.24%); Q= 1.982 MeV 

/?-; Q=11.508 MeV 

/?+ (99.89%), EC (0.11%); Q=2.754 MeV 
/3+; Q=5.145 MeV 
/?-; Q=0.156 MeV 



After the disappearance of one or two nucleons in the parent nuclide, one or 
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two holes appear in the nuclear shells; these holes will be filled in a subsequent 
nuclear de-excitation process, unless the nucleons reside on the outermost 
shells. If the initial excitation energy, E^xc, is higher than the binding energy 
Sn oi the least bound nucleon (A^ is p or n), the nucleus will be de-excited 
by particle emission (p, n, a, d, etc.); otherwise {Eexc < Sn) a 7 quantum 
will be emitted. In the following, we will take into consideration the and 
A^A^ decays from the last filled single-particle levels, when only 7 quanta could 
be emitted. Thus the daughter nucleus is exactly known. We also neglect by 
the prompt signal from gamma quanta emitted in initial deexcitation process 
(in particular, because of generally big half lives of daughter nuclei, and un- 
certainty in the number of emitted 7's, their energies, and rejection of high 
energy events by the muon veto), and search only for subsequent radioactive 
decay of created nuclei. 

The number of nucleons or nucleon pairs participating in the process can 
be calculated following refs. [13,15]. For example, after the decay of a neutron 
with binding energy E'^IA, Z), the excitation energy of the {A — 1,Z) daughter 
nucleus will be Eg^c = ^ni^^ ^) ~ Sn{A,Z). The condition to emit only 7 
quanta in the deexcitation process, Eexc < Sn{A — 1,Z), gives this restriction 
on the neutron binding energy: E^{A, Z) < Sn{A, Z) + SNiA — 1,Z). Similar 
equations can be written for p, pp, pn and nn decays (see ref. [15]). The values 
of the separation energies and S'jvtv were taken from ref. [34]. For single- 
particle energies of nucleons E^^^ on nuclear shells we used the continuum shell 
model calculations [35] for ^^C and ^^C, and the Hartree-Fock calculations with 
the Skyrme's interaction [36] for ^^O. 



3.2 Simulation of the response functions 

The expected response functions of the CTF detector and related efficiencies 
for the decay of unstable daughter nuclei were simulated with the EGS4 pack- 
age [37]. The number of initial electrons and 7 quanta emitted in the decay of 
the nucleus and their energies were generated according to the decay schemes 
[33] . The events were supposed to be uniformly distributed in the whole vol- 
ume of the hquid scintillator (and in a water layer close to the LS). The energy 
and spatial resolution of the detector [29] , light quenching factors for electrons 
and gammas [28], ^ veto and triggering efficiency were taken into account in 
the simulations. The calculated responses for the decay of ^^C and in the 
liquid scintillator (created after n and nn disappearance in ^^C, respectively) 
and ^^O in the water shield {nn decay in ^^O) are shown in fig. 2. In the last 
case only the water layer of 1 m thickness closest to the liquid scintillator was 
taken into consideration. Contributions from layers further out are negligible. 

The CTF background was simulated as well in order to check the understand- 
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Fig. 2. Energy distribution of the CTF-II installation collected during 29.1 d with 

all cuts. The expected response functions of the detector are also shown for ^^C 
(3.4- 10'^ decays in the liquid scintillator; corresponding mean lifetime for the n decay 
is Tn = 1.8-1025 y), (6.8- 10^ decays; r„n = 4.4-10^5 y), and ^'^O {lA-lO^ decays 
in 1 m thick water layer closest to the sphere with liquid scintillator; r„„ = 5.7 • 10^^ 

y)- 

ing of the detector. The simulated 1.46 MeV gamma peak together with 
experimental data is shown in the inset of fig. 1. The shift in the position of 
the full absorption peak is due to the ionization quenching effect (see section 
2.2). One can see the good agreement between the experimental and simulated 
peak position and resolution. 



3.3 Limits on probabilities of the N and NN disappearance 



The experimental data (see fig. 2, where the spectrum with all cuts is shown 
in detail) gives no strong evidence of the expected N and NN decay response 

functions, thus allowing only bounds to be set on the processes being searched 
for. In the present study, in order to extract the limits on the relevant mean 
lifetimes, we assumed conservatively that all events in the CTF experimental 
spectrum in the corresponding energy range AE are due to nucleon decays. 
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The mean lifetime limit was estimated using the formula 

'lim ^nucl ' -^obj ' ^/'S'lim -^nucl ' -^obj ' t/^liinn (1) 

where £ae is the detection efficiency in the AE energy window calculated 
in the full simulation of the relevant process, taking into account the radial 
cut efficiency en, and probability of identification by muon veto 77; Nnuci the 
number of parent nuclei; Nobj the number of objects {n, p or A'"A'" pairs) inside 
the parent nucleus whose decay will give the specific daughter nucleus; t the 
time of measurements; Sn^ the number of events (in the AE window) due to 
a particular effect which can be excluded with a given confidence level on the 
basis of experimental data; and = S\ua/sAE the corresponding number of 
decays in the liquid scintillator /water. 

The parameters in (1) were defined in the following way. The mass of scintil- 
lator was defined measuring the buoyancy with a precision of 5%. The total 
time of the data taking is 29.1 days and takes into account the dead time of 
the electronics of ~ 3%. The measurements were performed in runs of about 
24 hours. The internal source of '^'^K was used to check the energy calibra- 
tion stability. The position of the *^°K peak was defined for every run with 
statistical accuracy of ~ 20 keV. No systematical shift in energy scale with 
a time has been found. The same is valid for the energy resolution, which is 
directly connected to the energy scale. In addition, all the response functions 
for the searched decays are wider than detector's resolution, hence the result 
is insensitive to the observed energy scale uncertainties. 

The probabihty of detecting scintillation events by the muon veto -q — 1±0.2% 
has been defined from the experimental data for the energy = 1.9 MeV. 
This value is in agreement with the MC simulation. For 77(1. 9 MeV) = 1.2% 
the changes in the integral efficiency 1 — Tj{E) is negligible in the case of n and 
nn decays, for the p and pp decays the £ae will decrease by 7%. 

3.3.1 n and nn disappearance 

For the nn decay, there are four neutrons on the outermost 1^3/2 level of ^^C 
which gives the number of nn pairs Nobj = 2. The disappearance of the nn 
pair from this level will result in a ^'^C nucleus in the groimd state. Taking 
into account the statistical uncertainty in the number of experimental events 
in the 2.0-3.0 MeV energy window (276±17), we calculate Sn^a = 297 at 90% 
C.L. Accounting for the associated efficiency Eae — 0.44, the limiting value 
for decays is Dy,^ = 6.8 • 10^. With the values of N^^d = 1-9 ■ 10^9 for 
nuclei and t — 29.1 d, we obtain 

Tiim(rin,i2c) ^ 4 4 . IQ25 y ^j^ij c.L. 
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For the ^^O nucleus we will account for only one nn pair on the outermost lpi/2 
orbit {nn decay in deeper levels will result the ^^O nucleus being too excited). 
With the number of ^^O nuclei (in a 1 m thick water layer) N^uci = 9-8 • 10^^, 
Nohj — 1 and D^^ = 1.4 • 10^ for the energy region 2.2-2.6 MeV, the result is 

Tiim(rin/60) = 5.7 ■ 10^4 y with 90% C.L. 

With the assumption that the mean lifetime of the nn pair is the same in 
^^C and ^^O nuclei, and that nn decays in both of them contribute to the ex- 
perimental spectrum simultaneously, one can obtain a slightly more stringent 
limit for nn decay: 

Tiiminn inv) = 4.9 ■ lO^^ y with 90% C.L. 

For n decay in ^^C we used a similar approach, just demanding that the 
simulated response function for ^^C decay should be equal to the experimental 
spectrum in the energy region of 1.0-1.1 MeV (fig. 2). In this way the value 
Diim — 3.4 ■ 10^ was determined for the full number of ^^C decay events. 
Together with Ngbj — 4 (four neutrons on the lp3/2 level), it gives the following 
limit: 

Tiim(n,i2c) - 1.8 ■ 10^5 y with 90% C.L. 
3.3.2 p and pp disappearance 

As for the p and pp decays into invisible channels, the p disappearance in ^^C 

will result in ^^B nuclei, the f3~ decaying with high energy release (5=13.370 
MeV. The pp decays in ^"^C will produce ^^Be nuclei, also the f3~ decaying 
with (5=11.508 MeV (with probability of decay to the ground state of 57%). 

To estimate the Tiim for p and pp instabilities, we use the fact that no candidate 
scintillation events were observed in the CTF spectrum with energies higher 
than 4.5 MeV (fig. 2). 

High energy release in the liquid scintillator can activate the muon veto of the 
CTF, resulting in rejection of the event. We take into account such a suppres- 
sion of high energy tails in f3 decays of ^^B and ^^Be using the probability 1]{E) 
for identification of an event with energy E in the LS by the muon veto (fig. 3). 
The beta spectra of ^^B and ^^Be without and with suppression by the muon 
veto are also shown in fig. 3. The part of the ^^B beta spectrum with E > 4.5 
MeV reduced by a factor of l — r]{E) gives an integrated efficiency Eae — 0.39. 
With zero observed events (and with the assumption of zero expected back- 
ground), the limiting value for the number of events is Sum = 2.44 with 90% 
C.L. in accordance with the Feldman-Cousins procedure [38] recommended by 
the Particle Data Group [9]. Thus we arrive at the value Aim = 6.2 for ^^B 
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decay. Together with Nnud = 2.1 • 10^^ for parent ^^C nuclei and the number 
of A^obj — 4 for protons in orbit, we obtain 

riim(p/^C) = 1.1 • 10^6 y with 90% C.L. 
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Fig. 3. Probability of identification of an event with energy E in the scintillator 
by the muon veto (1). The (3 spectra of ^^B without (2) and with (3) suppression 
by the muon veto are also shown in arbitrary units. For signals with E < ~3 MeV 
(^"^C, ^^C and ^"^O decays), fi veto does not have a big effect on the overall efficiency. 
In the inset the experimental data taken with the radon source are presented. The 
tagged events of ^^^Bi-^^^Po at the energy £'=1.9 MeV are 'seen' by the muon veto 
system with an efficiency of = 0.01. 



In a similar way, for the pp decay in ^^C with the value Ngbj — 2 for the 
number of pp pairs and Eae — 0.36, the mean lifetime limit is 

Tiim(pp/^C) = 5.0 • 10^5 y with 90% C.L. 

All mean lifetime limits obtained here together with the numbers of parent 
nuclei, Nobj and the numbers of decay events are summarized in table 3. 
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Table 3 

Mean lifetime limits, th^, (at 90% C.L.) for N and NN decays in the CTF. Nnud 
is the number of parent nuclei; Ngbj the number of objects (n, p and NN pairs) per 
parent nucleus; -Dum the excluded number of decay events. 



Decay Nnud -/Vobj -Diim Tum, y 



p 




2.1 


10^7 


4 


6.2 




1.1 


1026 


n 


12p ,llp 


1.9 


1029 


4 


3.4 


10^ 


1.8 


1025 


nn 




1.9 


1029 


2 


6.8 


102 


4.4 


1025 






9.8 


1029 a 


1 


1.4 


10^ 


5.7 


1024 


pp 


ifC^ilBe 


2.1 


1027 


2 


6.7 




5.0 


1025 



" In 1 m thick layer of water closest to the CTF liquid scintillator 
4 Conclusions 



Using the unique features of the BOREXINO Counting Test Facility - the 
extremely low background, the large scintillator mass of 4.2 ton and the low 
energy threshold - new limits on N and NN decays into invisible channels 
(disappearance, decays to neutrinos, majorons, etc.) have been set: 

r{n inv) > 1.8 • 102^ y, r(p inv) > 1.1 ■ 102*^ y, 

T{nn inv) > 4.9 • 102^ y and r{pp inv) > 5.0 ■ 102^ y with 90% C.L. 

Comparing these values with the data of table 1, one can see that the estab- 
lished limits for nn and pp decays are the best up-to-date limits set by any 
method, including radiochemical and geochemical experiments. 

These limits are obtained in a very conservative assumption that all events 
in the corrresponding energy region are due to the nucleon decays. The data 
from the full scale Borexino detector will improve the presented limits by at 
least two orders of magnitude. 
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